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Abstract 

Universal aspects of the thermodynamic Casimir effect in wetting films of 3 He- 4 He mixtures 
near their bulk tricritical point are studied within suitable models serving as representatives of 
the corresponding universality class. The effective forces between the boundaries of such films 
arising from the confinement are calculated along isotherms at several fixed concentrations of 3 He. 
Nonsymmetric boundary conditions impose nontrivial concentration profiles leading to repulsive 
Casimir forces which exhibit a rich behavior of the crossover between the tricritical point and 
the line of critical points. The theoretical results agree with published experimental data and 
emphasize the importance of logarithmic corrections. 

PACS numbers: 05.70.Jk, 64.60.Fr, 64.60.Kw, 67.40.Kh, 68.15,+c 
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Finite-size contributions to the free energy of a fluid confined between two surfaces at 
a distance L give rise to an effective force between them. Theory predicts that at the 
bulk critical point T c of such a system this force becomes long-ranged as a result of critical 
fluctuations of the corresponding ordering degrees of freedom. This is analogous to the well- 
known Casimir effect in electromagnetism. This so-called critical Casimir force /c_P er un it 
area and in units of ksT c can be expressed in terms of universal scaling functions jl|. 

Only recently, sophisticated wetting experiments have provided detailed quantitative data 



for critical Casimir forces in various systems . In the case of He wetting films 

near the superfluid transition, these experimental studies support quantitatively theoretical 
predictions for fc (T > T c ) corresponding to the universality class of the XY model 
For the case of 3 He- 4 He films near the bulk tricritical point some theoretical predictions are 

Q 

available [2fl, but those do not apply for the boundary conditions relevant for recent wetting 
experiments performed in these systems "\ \^ . However, the shape of the scaling function of 
the Casimir force depends sensitively on the type of boundary conditions (BC) and thus on 
the surface universality classes to which the confining surfaces belong 8]. The experiments of 
Re, Q rep0rt a ^ h around the tricritica! pent which SUgg ests nonsymmetric BC for 
the superfluid order parameter (SOP). This is opposite to the case of pure 4 He wetting films 
near the A-point where fc was found to be attractive . For the latter system the BC seem 
to be very well approximated by symmetric Dirichlet-Dirichlet BC (O, O) forming the so- 
called ordinary (O) surface universality class because the quantum-mechanical wave function 
describing the superfluid state vanishes at both interfaces 

ELQ 

. The type of BC for 3 He- 4 He 

wetting films is not clear from the outset because a 4 He-rich layer forms near the substrate- 
fluid interface, which may become superfluid already above the bulk A-line Q| whereas 3 He 
has a preference for the fluid-vapor interface. Thus the two interfaces impose a nontrivial 
concentration profile which in turn couples to the SOP. This leads to the hypothesis that the 
concentration profile induces effectively nonsymmetric (O, +) BC for the SOP, i.e., Dirichlet 
boundary conditions at the fluid-vapor interface and symmetry-breaking (+) BC at the 
substrate-fluid interface (also known as the so-called extraordinary or normal universality 
class j8j). For the present tricritical behavior the upper critical dimension d* equals 3. 
In this case theory predicts that for three-dimensional systems the asymptotic tricritical 
thermodynamic functions exhibit power laws with critical exponents taking their classical 
values. However, logarithmic corrections to the mean-field (MF) behavior are expected 
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under experimental conditions [lo| . 

Here we consider two complementary approaches. Field-theoretical methods and 
renormalization-group (RG) analyses are used to derive universal properties of the Casimir 
force at the tricritical point and the form of logarithmic corrections. However, these methods 
do not lend themselves for systematic studies of fc along all thermodynamic paths followed 
in the aforementioned experiments. In order to be able to interpret the rich variation of fc 
extracted from the capacity measurements in Ref. k|, to understand the emergence of the 
actual BC and, moreover, to predict the behavior of fc in the crossover region between the 
tricritical and the critical points, we employ the vectoralized Blume-Emery-Griffiths model 
(VBEG) |l2j as a representative of the same universality class as the actual physical system. 
This lattice model is extended to the film geometry and treated within mean field theory 
(MFT). 

First we derive the leading asymptotic behavior of fc at tricriticality for (O, +) BC. 
To this end we consider the standard Ginzburg-Landau (GL) Hamiltonian for an 0(n)- 
symmetric tricritical system (T = T t ) in a film geometry: 

«[*] = / d^x £ dz + (i) 

where L is the film thickness, $ is the n-component OP and z is the distance between the 

confining surfaces; u is a bare coupling constant. In a film geometry fc = —(df ex /dL) = 

(T zz ) is given by the stress tensor component T zz [1], where f ex (L) = (/ — fb)L, f is the 

total free energy per unit area and per ksTt and /& is the bulk contribution. The stress 

tensor is given by [l] %j = ■ dj& — 5ijC - (d — 2)/(4(d— l))(didj - (%V 2 )<fr 2 , where £ is 

the integrand of (0). We take $ = (m(z), 0, . . . , 0). Determination of the tricritical Casimir 

force starts from the Euler-Lagrange equation for the OP profile: m"(z) = (u/120)m 5 (z) 

with (O, +) BC, i.e., m(0) = and m(L) = +oo. In this case the spatially constant (% z ) 

can be expressed as (l/2)(m'(0)) 2 . With the scaling ansatz m(z) = (u / '360) -1 1 ' 4 L'^^z/L) 

and (T zz ) = (90/u) 1//2 L _3 O, and after integrating directly the first integral of the Euler- 

Lagrange equation one obtains O 1 / 3 = / dp/ yl + p 6 ~ 1.40218 by implementing the 

Jo 

above BC. Eventually, in units of ksT t the MFT result for the tricritical Casimir force f c in 
the case of (O, +) BC is f c = 2.7568(4) (90/n) 1/2 L~ 3 . Note that within MFT the parameter 
u > remains undetermined. Its value follows from using standard RG arguments. In 
d = 3 — e the above MFT result yields the leading contribution in an e-expansion. After 



removing the uv singularity via renormalization the asymptotic scaling behaviour of /£, 
follows from substituting u by the appropriate fixed-point value m* oc e. At d — d*, and 
under spatial rescaling by a factor £*u flows to its RG fixed point value u* = according 
to u(£) = (240tt 2 )/ ((3n + 22)\\n£\) [ll|. With the rescaling factor i = l Q /L, where l is a 
microscopic length scale of the order of a few A, this yields a logarithmic correction to the 
power law L-dependence of the tricritical Casimir force: 



f c ~ 0.54(3n + 22) 1/2 (ln L/l ) 1/2 L~ 3 . (2) 

xbc 

L 



d dominant, 
for their 



Gaussian fluctuations give contributions of at least 0(u°) which are therefore su 
We compare Eq. (J2J) for n = 2 with the data obtained by Garcia and Chan 
experimental value L/l pa 520 A/1.3 A. This gives "&t = fcL 3 pa 6.96 in a good agreement 
with d e t xp = 8.4 ± 1.7, which suggests that this experiment maybe the first to have verified 
implicitly the existence of logarithmic corrections near the tricritical point. However, in 
order to extract the actual value of the universal Casimir amplitude (i.e., the numerical 
prefactor in Eq. (J2J) the experimental data require a reanalysis based on a functional form 
given by Eq. (J2J). 

Now we turn to the VBEG model and consider a d — 3 simple cubic lattice consisting of L 
parallel lattice layers at spacing a. Each layer has A sites, labeled i,j, . . . and associated with 
an occupation variable ti — 0,1 and a phase Q{ (0 < 0j < 2tt) which mimics the phase of the 
4 He wave function. A 3 He ( 4 He) atom at site % corresponds to ti = 0(1). The Hamiltonian 
is given by 

<ij> <ij> i 

where the first two sums run over nearest-neighbor pairs and the last one is over all lattice 
sites. The field Aj is related to the chemical potentials of the two components of the mixture. 
Aj = Ai and Aj = A2 on the left and right surface layer, respectively, and Aj = A otherwise. 
The differences Aj — A, i = 1,2, are a measure of the relative preferences of 4 He atoms for 
the two surfaces such that Ai < A corresponds to the preference of 4 He atoms for the solid 
substrate. At the opposite surface with the vapor we choose A 2 = A t , the tricritical bulk 
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value. Near the tricritical point this choice is consistent with the assumption made in Ref. 
for the concentration profile across the wetting film, whereby at the interface with the vapor 
the 3 He mole fraction takes the bulk value. Bulk properties of this model were studied 



within MFT and by Monte Carlo simulations in d = 3 in 



12j. The resulting bulk phase 



diagram resembles that observed experimentally for 3 He- 4 He mixtures (see Fig. [TJ. For the 
film geometry we have solved this model within MFT. This yields a set of self-consistent 
equations for the OP in the /th layer, i.e., the concentration Qi =< t\ >= 1 — X(l) of 
4 He, and the SOP = (m^^m^) with m z = (ticosOi) and m| = (i^sin^) where ti 
and 9i denote the occupation number and the phase in the /th layer, respectively. In the 



absence of helicity one has = (m[ 1} ,0) = (m ; ,0), Qi = I (Jbi)/ ( e -^ a;+Ai/T + I (Jbi) 



and A; = A for I ^ 1, L; K — qK/T and J = qJ/T, where q = q± + 2q' is the coordination 
number of the lattice. q± is the in-layer coordination number while each site (but not in the 
first and last layers) is connected to q' atoms in each adjacent layer. We have introduced 
b\ = mi_i+qi_mi+mi + i for I 1, L, b\ = q±m 1 +m 2 , and b L = m L _i+q±m L and analogously 
ai = Qi-i + q±Qi + Qi + i for I ^ 1, L, ai = q±Qi + Q 2 , and a L = Q L -i + q±QL- The coupled 
equations for Q\ and m\ are solved numerically; the acceptable solution minimizes the free 
energy. First, we have analyzed the semi-infinite system. Close to the A-line we observe 
a higher 4 He concentration near the left surface, which induces a local superfluid ordering 
(see Fig. |2J). By varying T and A one obtains a whole line of continuous surface transitions 
corresponding to the onset of the formation of a superfluid film near the wall; it meets the 
A-line at the special transition point whose position depends on the value of Ai (see Fig. P). 
These findings are in agreement with the results of a Migdal-Kadanoff analysis |l3j |. In the 
film geometry the Casimir force is obtained by taking a finite difference after calculating f ex 
for Lq and Lq + 1. Figure El summarizes our result for a film of width L = 20, K/J = 0.5, 
Ai/J = —3, and A 2 = A t /J ~ 0.61. fc is calculated along the thermodynamic paths 
indicated in Fig. ^ Below T t , fc is calculated along the coexistence line, infinitesimally 
on the superfluid branch of bulk coexistence. Our results are presented in terms of the 
scaling function = L d fc as a function of the scaling variable y = tL 1 ^ = (L^o/0 1 ^ > 
where t = (T — T t )/T t . £ is the amplitude of the correlation length £ and v = 1. The 
surface transition does not leave a visible trace in the behavior of fc- For 3 He concentration 
X < X t , upon crossing the A-line there is a steep variation associated with a break in slope, 

[iving rise to the formation of shoulders which are similar to those observed experimentally 
|. For X > X t , when T reaches the phase separation temperature fc coincides with the 
curve common to all values of X. This occurs with a discontinuous first derivative. The 
aforementioned common curve exhibits a pronounced maximum below T t at y ~ —0.74 and 





Iq(z) and h(z) are modified Bessel functions 
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gradually decreases to zero for y — > — oo. Below T f , both the concentration and the SOP 
profiles corresponding to this common curve display an interface-like structure separating 
two domains of the coexisting bulk phases (see t = —0.0625 in Fig. |2J). Features of fc in 
this 'soft mode' phase can be attributed to purely interfacial effects, similarly to Ising-like 
films with asymmetric BC [14]. Beyond MFT a positive sign of the force can be regarded 
as a consequence of entropic repulsion [15] . The maximum of fc is expected to occur at the 
temperature T for which the interfacial width ~ £5 ~ L, i.e., y ~ — 1 [lJJ which checks with 
Fig. |3] Q]. For X < X t — 0.05 we observe a crossover to the critical superfluid behavior 
of pure 4 He and a gradual formation of a second, less pronounced local maximum located 
slightly below the A-line. This local maximum decreases upon departure from X t and finally 
disappears above the special transition S. This is expected, because above S the BC turn 
into the type (O, O) for which fc vanishes within MFT. For lower T, fc increases steeply 
upon approaching bulk coexistence revealing that interfacial effects associated with the 'soft 
mode' lead to a much stronger Casimir effect than near the critical A-line. The qualitative 
features of $ extracted from the experimental data for X ~ X t (see Fig. 5 in Ref. 
are very well captured by the present lattice model. Discrepancies can be attributed to 
fluctuation effects neglected in the present MFT VBEG approach: (i) The discontinuities of 
slopes as obtained within MFT upon crossing the A-line are expected to be smeared out by 
fluctuations, (ii) The experimental scaling function d does not vanish at low temperatures, 
which may be due to Goldstone modes in the superfluid phase. However, the possibility that 
this behavior is an artifact of an extreme change in the dielectric constant of the film cannot 



be excluded 



171 ] . In the crossover regime to the critical behavior only few experimental 



data for the thickness of the wetting film are published. However, again the variation of 
film thicknesses agrees with our findings. In particular, one observes a rapid thickening of 
the films upon approaching the coexistence line; for some values of X a small maximum 
located near the A-line is also visible (compare Fig. EJ). A quantitative comparison is not 
possible because, for our choice of surface terms in the Hamiltonian, the fixed-point BC 
(O, +) cannot be reached within the VBEG model. In order to be able to extract universal 
properties - which requires to reach fixed-point BC - it would be necessary to introduce a 
surface field which couples to the SOP so that the BC (+) can be realized. Also MFT is not 
sufficient, even in d — d*. A naive correction of •& by the logarithmic factor derived within 
the GL model will not give the proper universal behavior. Instead, renormalization group 
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schemes or Monte Carlo simulations have to be employed. Nonetheless, our MFT results for 
X = X t , if matched at the tricritical point y = and after adjusting the amplitude £o so 
that the positions of the maximum of the scaling function are the same (i.e., £q /a ~ 0.065) 
reproduce very well the experimental curve (see Fig.0J), especially near the maximum where 
we expect interfacial effects to be dominant. This is consistent because the 'soft mode' 
phase does not depend on the details of the surface fields. Notice, that the experimental 
data nominally for X = X t more closely match the theory for X = X t — 0.01. This raises 
the question as to whether the 3 He concentration in the film is shifted relative to the bulk 
one. 

A.M. benefited from discussions with R. Garcia, M. Krech and S. Kondrat. This work 
was partially funded by KBN grant No.4 T09A 066 22. 
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FIG. 1: Phase diagram for the VBEG model obtained whithin MFT for K/J = 0.5 and A\/J = —3 
exhibiting the bulk A-line T S (X) of continuous superfluid transitions (long-dashed line), the phase 
separation curves (solid lines), the tricritical point A = (T/T s (0) = 2/3, X = 1/3), and the surface 
transition line (short-dashed line) which merges with the bulk A-line at the special transition point 
S = (T/T s (0) ~ 0.759, X ~ 0.241). Vertical lines represent thermodynamic paths along which the 
Casimir force has been calculated (see Fig. EJ). ♦,• (A), ▲: state points considered in Fig. |2 
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FIG. 2: (a) 3 He concentration profile X(l) = 1 — Qi and (b) SOP mi profile for a VBEG film of 
width L = 60 for K = 0.5 J, A±/J = —3, and A2/J = At/ J ~ 0.61 corresponding to the state 
points ♦, •, and ▲ indicated in Fig. ^ 
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y = tL 

FIG. 3: Color online. Dimensionless scaling function i}(y = tL) = fcL^ for the VBEG model 
calculated within MF theory along the paths of fixed concentration of 3 He shown in Fig. 1. Dots 
indicate the corresponding onset temperature T s (y) of superfluidity at the A-line. 
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y = tL 



FIG. 4: The adjusted scaling function (see main text) for the VBEG model within MFT 



compared with the corresponding experimental curve 
concentration Xt ~ 0.672 of 3 He. 



J] obtained along the path of fixed tricritical 
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